Detection of low Eu abundances in extremely metal-poor stars 
and the origin of r-process elements 1 
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ABSTRACT 

We report abundance analyses of three extremely metal-poor stars with [Fe/H] < — 3, using the 
Subaru High Dispersion Spectrograph (HDS). All are found to have sub-solar values of [Eu/Fe]. 
Comparison with our chemical evolution model of the Galactic halo implies the dominant source 
of Eu to be the low-mass end of the supernova mass range. Future studies of stars with low Eu 
abundances will be important to determine the r-process site. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: Popula- 
tion II — supernovae: general — Galaxy: evolution — Galaxy: halo 



1. Introduction 

The origin of the rapid neutron capture (re- 
process) elements is still uncertain. Although a 
few scenarios such as neutrino winds (Woosley et 
al. 1994) in core-collapse supernovae (SNe), the 
collapse of O-Ne-Mg cores resulting from 8 — 10M Q 
stars (Wanajo et al. 2003), and neutron star merg- 
ers (Freiburghaus et al. 1999) show some promise, 
no consensus has been achieved. 

Analyses have shown a large scatter in the 
abundances of neutron-capture elements such as 
Sr, Ba, and Eu among metal-poor stars in the 
Galactic halo (e.g., McWilliam et al. 1995; Ryan, 
Norris, & Beers 1996). If the scatter accurately 
portrays the state of the interstellar medium 
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(ISM) from which they formed, large dispersions 
in excess of observational errors indicate that the 
ISM was not fully mixed at the early epoch of 
Galactic history and that metal-poor stars con- 
tain products from only one or a few SNe (Ryan, 
Norris, & Bessell 1991). The scatter possibly re- 
flects variation in the yields of SNe from different 
mass progenitors; the huge dispersions in neutron- 
capture elements imply that r-process yields are 
highly dependent on the masses of SN progenitors. 

Eu traces the r-process. The chemical evolu- 
tion model of Ishimaru & Wanajo (1999, hereafter 
IW99) explains the large dispersion of Eu/Fe ob- 
served in very metal-poor stars if the r-process 
elements originate from a limited mass range of 
progenitor stars, such as the SNe from stars of 
8 — IOMq or alternatively > 30M Q (see also 
Travaglio et al. 1999). Tsujimoto, Shigeyama, 
& Yoshii (2000) have concluded by a similar ap- 
proach that the SNe from ss 2OM are the domi- 
nant source of the r-process elements. A clear dif- 
ference among these cases is expected in the stellar 
distribution of [Eu/Fe] 6 at [Fe/H]< -3 (IW99). 
However, a shortage of observational data in this 
metallicity range has made it difficult to distin- 
guish between the proposed r-process sites. 



"^[Xi/Xj] = \og(Ni/Nj) - log(Ni/Nj) e , where N t indicates 
abundance of i-th element X;. 
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In this Letter, we report on three extremely 
metal-poor stars which we show to have very low 
Eu abundances (§2). These data are compared 
with our chemical evolution models to distinguish 
between the proposed r-process sites (§3). 

2. Observations and Analysis 

We selected three very metal-poor ([Fe/H] 
< -3) giants, HD 4306, CS 22878-101, and 
CS 22950-046, which were known from previous 
studies (McWilliam et al. 1995; McWilliam 1998) 
to have [Ba/Fe] ~ — 1, typical for their metallici- 
ties. Observations were made with the High Dis- 
persion Spectrograph (HDS; Noguchi et al. 2002) 
of the 8.2m Subaru Telescope in 2001 July, at a 
resolving power R = 50, 000. The total exposure 
time was 50, 240, and 360 minutes for HD 4306, 
CS 22878-101, and CS 22950-046, respectively, 
resulting in S/N ratios per 0.012 A pixel at 4100 
A of 260, 110, and 80. Data reduction was per- 
formed in the standard way with IRAF 7 . 

Equivalent widths were measured by fitting 
Gaussian profiles to the absorption lines of species 
listed in Table l 8 , except for Eu. An LTE anal- 
ysis using model atmospheres of Kurucz (1993), 
based on equivalent widths measured above, was 
performed for species other than Eu. Results are 
given in Table 1. For Ba we included the hyperfinc 
splitting and isotope shifts (McWilliam 1998), as- 
suming isotope ratios of the r-process component 
in solar system material (Arlandini et al. 1999). 

The atmospheric parameters (effective temper- 
ature: T c ff (K), gravity: g (cm s~ 2 ), micro- 
turbulent velocity: u m i c (km s _1 ), and metallic- 
ity that is assumed to be [Fe/H]) of CS 22878- 
101 and CS 22950-046 were adopted from Co- 
hen et al. (2002) and Carretta et al. (2002), i.e., 
T cff /log 5 /[Fe/H]/zw - 4775/1.3/-3.1/2.0, and 
4730/1.3/-3.3/2.0, respectively. We found that 
these gravities satisfy the ionization balance be- 
tween Fe I and Fe II, and the microturbulence 
velocities lead to no dependence of the derived 
abundance on equivalent widths of Fe I lines. 

7 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 

8 The line list and measurements arc available at 
http://optik2.mtk.nao.ac.jp/~waoki/EW/ 



The derived Fe abundances also agree very well 
with the assumed metallicities. So, we adopted 
these atmospheric parameters without modifica- 
tion. For HD 4306, we adopted T cff = 5000 K 
based on V - K = 2.15 (Norris, Bessell, & Pick- 
les 1988; and 2MASS data: Skrutskic et al. 1997) 
and the temperature scale of Houdashelt, Bell & 
Sweigart (2000). We derived g and w m i c from 
the analysis by requiring ionization balance and 
no dependence of abundance on line strength. 
The resulting parameters, T e ff/logg/[Fe/H]/i; m i c 
= 5000/2. 0/-3. 0/1.85, are similar to those of 
other works (McWilliam et al. 1995, Honda et al., 
in preparation). 

For the three Eu II lines, at 3819, 4129, and 
4205 A, a spectrum synthesis technique was ap- 
plied. The line list includes hypcrfme splitting and 
isotope shifts from Lawler et al. (2001). Contam- 
ination by other species was included using the 
comprehensive line list of Kurucz (1995). Figure 1 
compares the observed and synthetic spectra for 
the 3819 A line, which is the strongest of the three. 
Eu is clearly detected in the spectra of HD 4306 
and CS 22878-101 but not in CS 22950-046. We 
derived consistent results from the 4129 A line, but 
the measurement of 4205 A is uncertain because 
of severe contamination by V II at 4205.09 A. 
We adopted the average of the Eu abundances de- 
rived from the 3819 and 4129 A lines. In Table 2, 
we give the equivalent widths of the Eu II lines 
estimated from the spectrum synthesis using the 
Eu abundances derived from individual lines. For 
CS 22950-046, we estimated the upper limit of its 
equivalent widths from the 3 a depth of the ran- 
dom noise level and typical line widths for this 
object, and determined the upper limit of the Eu 
abundance, given in Table 1. The continuum fit 
is constrained over a wider spectral range that is 
shown in Figure 1; a systematic error of 1% would 
change the [Eu/H] limit from the 3819 A line by 
0.1. 

The random error for each species is estimated 
from the standard deviation of the abundances de- 
rived from individual lines. The values are some- 
times unrealistically small when the number of 
lines is small, so we adopted the larger of the value 
for the species and that for Fe I, 0.08-0.15 dex 
depending on the species and the quality of the 
spectra. We estimated the errors in the Eu abun- 
dances by eye from fitting synthetic spectra to 
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the observed ones, adopting 0.1 dex for HD 4306 
and 0.2 dex for CS 22878-101. The errors in 
the abundance determinations due to uncertain- 
ties of the atmospheric parameters were evaluated 
for AT eff = 100 K, Alogg = 0.3, Av mic = 0.3 km 
s-\ and A [M/H] = 0.3 for HD 4306, and added 
in quadrature with the random error estimated 
above, to give the total uncertainty in Table 1. 
The larger errors in [Sr/Fe] reflect the sensitiv- 
ity to microturbulence of the two resonance lines. 
The [Eu/Fe] values are affected most significantly 
by the uncertainty in the surface gravity, by 0.12 
for Alogg = 0.3. 

Comparisons between our results for the two CS 
objects and those by Carretta et al. (2002) show 
good agreement (within 0.17 dex) for the species 
in Table 1 other than Eu. An exception is Ba in 
CS 22878-101: our abundance is 0.35 dex higher 
than that by Carretta et al. (2002). The reason for 
this discrepancy is unclear. The equivalent widths 
of the Ba II 4554 A line agree well between the two 
studies, and the atmospheric parameters are the 
same. The difference in the assumed isotope ratios 
of Ba in the abundance analyses do not explain 
this discrepancy. 

3. Comparison with the Galactic chemical 
evolution models 

Assuming that star formation is induced by 
a single supernova explosion, IW99 have con- 
structed an inhomogcncous chemical evolution 
model of the Galactic halo. The chemical com- 
position of a newly formed star is a mixture of 
the products from each SN and the ISM swept up 
by the expanding ejecta. This model, along with 
some improvements (see Ishimaru, Prantzos, & 
Wanajo 2003), leads to a large scatter of relative 
abundances for elements that have strong depen- 
dencies on the progenitor masses. The scatters of 
relative abundances in a and iron-peak elements, 
on the other hand, are predicted to be small (as 
is observed) owing to the weak dependencies of 
their yield ratios on the progenitor masses. We 
use this model to derive the Galactic evolution of 
the r-process tracer Eu. We consider three cases, 
where the r-process elements originate from SNe 
of (a) 8 - 1OM , (b) 20 - 25M , and (c) > 30M Q 
stars. 

Yields for Type II and Type la SNe are taken 



from Nomoto et al. (1997a, b). The 8-10M Q stars 
are assumed to produce no iron, since their contri- 
bution to the enrichment of iron-peak elements in 
the Galaxy is negligible (Wanajo et al. 2003). The 
mass of Eu produced per SN event is assumed to 
be constant over the mass range of each case; this 
simplification will not affect the result significantly 
because of the narrow mass range of Eu-generating 
SNe assumed in each case. The requirement that 
the model reproduces the solar values [Eu/Fe] = 
[Fe/H] = implies ejected Eu masses 3.1 x 10~ 7 , 
1.1 x 10~ 6 , and 7.8 x 1O~ 7 M for cases (a)-(c), re- 
spectively. These are reasonable amounts from a 
nucleosynthetic point of view (Wanajo et al. 2001, 
2003). ' 

Figure 2 compares the results of these cases 
with the observational data from §2 (large dou- 
ble circles) and from sources cited in the caption 
(small circles). The predicted log number den- 
sity of stars per unit area on the [Fc/H]- [Eu/Fe] 
plane is also shown as a color image. The thick 
solid lines indicate the average values of [Eu/Fe] 
as a function of metallicity, calculated from the 
models. The solid and thin solid lines indicate, 
respectively, the 50% and 90% confidence lines. 

Figure 2 shows that observable differences be- 
tween the cases appear at [Fe/H]< —3. In case (c), 
most of the stars are expected to have [Eu/Fe] > 0, 
owing to Eu production solely by massive, short- 
lived stars. In cases (a) and (b), significant num- 
bers of stars having [Eu/Fe] < are predicted at 
[Fe/H] < —3 owing to the delayed production of 
Eu by lower mass SN progenitors. This is out- 
standing in case (a), which has the largest delay. 

Most previous observational data, many of 
them upper limits, distribute between the 90% 
confidence lines for all cases, which has made it 
difficult to determine the mass range of the r- 
process site. Our newly obtained data add the 
lowest detections of Eu, at [Fe/H] < —3, and help 
distinguish between the three cases. The best 
agreement can be seen in case (a), in which the 
three stars, and most other stars from previous ob- 
servations are located between the 50% confidence 
lines at [Fe/H] < —3. In case (b), all three stars 
are located below the average (thick solid) line, 
although they remain above the 90% confidence 
line. In case (c), these stars are located outside 
the 90% confidence region. We suggest, therefore, 
that case (a) is most likely to be the r-process 
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site, i.e. SNe from low-mass progenitors such as 
8 — 10M Q stars. More detections of Eu at low 
metallicity will be needed to clearly distinguish 
between cases (a) and (b). Note that this result 
is not due to a selection bias, since the [Ba/Fe] 
values (~ — 1) are typical of extremely metal-poor 
stars with [Fc/H] 3 (Ryan et al. 1996). 

Our analysis gives [Ba/Eu] values consistent 
with the solar r-process (Arlandini et al. 1999) 
when estimated errors are included (see Table 3). 
Hence our result may hold for heavy r-process el- 
ements with Z > 56, not just Z ~ 63. The values 
of [Sr/Ba], however, are significantly higher than 
for the solar r-process, implying that these three 
stars exhibit light r-process elements (Z < 56) 
produced in more massive SNe (> 10M Q ). 

The discussion above suggests that the produc- 
tion of the r-process elements is associated with a 
small fraction of SNe near the low-mass end of the 
range. Neutrino winds in the explosions of mas- 
sive stars may face difficulties in being a dominant 
source of the r-process elements. Wanajo et al. 
(2001) have demonstrated that an r-process in the 
neutrino winds proceeds from only very massive 
proto-neutron stars, which might result from mas- 
sive progenitors such as > 20 — 30M Q stars, which 
is similar to case (b). Hypernovae (> 20 — 25M , 
Maeda & Nomoto 2003) or pair-instability super- 
novae (140 - 26OM , Heger & Woosley 2002), re- 
sulting from stars near the high mass-end of the 
SN progenitors, similar to case (c), are clearly ex- 
cluded as the major r-process site. 

We suggest, therefore, that the dominant r- 
process site is SN explosions of collapsing O-Ne- 
Mg cores from 8 — 10M Q stars (Wheeler, Cowan, & 
Hillebrandt 1998). Recently, Wanajo et al. (2003) 
have demonstrated that the prompt explosion of 
the collapsing O-Ne-Mg core from a 9M Q star 
reproduces the solar r-process pattern for nuclei 
with A > 130, and is characterized by a lack of 
a-elements and only a small amount of iron-peak 
elements. This clearly differs from more massive 
SNe with iron cores (> 1OM ) that eject both 
these elements, and is consistent with the fact 
that the abundances of heavy r-process elements 
in stars with [Fe/H] ~ — 3 are not related with 
those of iron-peak elements or of elements with 
lower atomic numbers (Qian & Wasserburg 2003). 

This study shows the importance of detecting 
Eu in extremely metal-poor stars to explore the 



origin of r-process elements. Further observations 
are needed to confirm the origin. 
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Table 1 
ELEMENTAL ABUNDANCES 







HD 4306 






CS 22878-101 






CS 22950-046 






[X/Fc] a 


n 


a b 


loge 


[X/Fe] a 


n 


a b 


loge 


[X/Fe] a 


n 


a b 


loge 


Fel 


-2.76 


168 


0.13 


+4.74 


-3.14 


161 


0.14 


+4.36 


-3.34 


141 


0.18 


+4.16 


Fe II 


-2.79 


19 


0.13 


+4.71 


-3.07 


20 


0.14 


+4.43 


-3.34 


17 


0.18 


+4.16 


Sr II 


+0.28 


2 


0.27 


+0.42 


-0.15 


2 


0.28 


-0.33 


-0.18 


2 


0.30 


-0.60 


Ba II 


-1.09 


2 


0.14 


-1.65 


-0.73 


2 


0.16 


-1.61 


-1.31 


2 


0.19 


-2.43 


Eu II 


-0.57 


2 


0.16 


-2.80 


-0.30 


2 


0.24 


-2.85 


< -0.2 


1 




< -3.C 



a [Fe/H] for Fe I and Fe II. 
b Uncertainty in [Fe/H] or [X/Fe] values. 
c loge(X) = \og w {N x /Xy) + 12.0 for elements X. 



Table 2 

EQUIVALENT WIDTH OF EU LINES 



Wavelength 




Equvalcnt Width (mA) 


A 


HD 4306 


CS 22878-101 CS 22950-046 


3819.7 


3.8 


7.4 <6.8 


4129.7 


2.2 


5.3 


4205.0 


2.5 


5.0 <3.4 



Table 3 

Abundance rations among Sr, Ba, and Eu 





HD 4306 


CS 22878-101 


CS 22950-046 


solar''' 


[Sr/Ba] 


1.37 + 0.30 


0.58 + 0.32 


1.13 + 0.36 


-0.10 


[Ba/Eu] 


-0.52 + 0.21 


-0.43 + 0.29 


> -1.11 


-0.69 



'''solar system r-process abundances by Arlandini ct al. (1999) 
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0.9 


-CS22878- 


101 [Eu/Fe] = 


-0.35- 
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-CS22950- 
i 


046 [Eu/Fe]< 
i 


-0.2 " 
i 



3819.2 3819.6 

Wavelength (A) 



Fig. 1. — Comparison of the observed spectra 
(dots) and synthetic ones (lines) near the Eu II 
3819.7 A line. For HD 4306 and CS 22878-101, 
three synthetic spectra differ in step of A[Eu/Fe] 
= 0.2 dex are shown by solid lines. Central line is 
calculated for the Eu abundance presented in each 
panel. The solid line in the panel for CS 22950- 
046 indicates the synthetic spectrum for the upper 
limit of Eu abundance presented in the panel (see 
text). Dotted lines are the spectra calculated as- 
suming no Eu contribution. 




[Fe/H] 



Fig. 2. — Comparison of the observed data with 
the model predictions. The r-process site is as- 
sumed to be SNe of (a) 8 - 1OM , (b) 20 - 25M , 
and (c) > 30Mq stars. The predicted number 
density of stars per unit area is color-coded. The 
average stellar abundance distributions are inde- 
cated by thick-solid lines with the 50% and 90% 
confidence intervals (solid and thin-solid lines, re- 
spectively). The average abundances of the ISM 
are denoted by the thick-dashed lines. The current 
observational data are given by large double cir- 
cles, with others (small circles) taken from Grat- 
ton & Sneden (1994); McWilliam et al. (1995); 
McWilliam (1998); Woolf, Tomkin, & Lambert 

(1995) ; Ryan, Norris, & Beers (1996); Shetrone 

(1996) ; Sneden et al. (1996); Westin et al. (1998); 
Burris et al. (2000); Fulbright (2000); Norris, 
Ryan, k Beers (2001); Johnson (2002); Johnson 
& Bolte (2002); Frangois et al. (2003); Honda (pri- 
vate communication). 
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